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Abstract. Several new formulas are developed that enable the evaluation of a family of definite integrals 
containing the product of two Whittaker W K , M (a;)-functions. The integration is performed with respect to 
the second index fi, and the first index k is permitted to have any complex value, within certain restrictions 
required for convergence. The method utilizes complex contour integration along with various symmetry 
relations satisfied by the Whittaker functions. The new results derived in this paper are complementary to 
the previously known integrals of products of Whittaker functions, which generally treat integration with 
respect to either the first index k or the primary argument x. A physical application involving radiative 
transport is discussed. 



I. INTRODUCTION 

Many problems in mathematical physics involve differential equations with solutions that can be expressed 
in terms of Whittaker 's functions W K:fJ ,(x) and M K ^{x). Examples of the diverse applications include studies 
of the spectral evolution resulting from the Compton scattering of radiation by hot electrons, 1 ' 2,3 modeling of 
the structure of the hydrogen atom, 4 analysis of the Schrodinger equation, 5 studies of the Coulomb Green's 
function, 6 and analysis of fluctuations in financial markets. 7 

In a number of applications, it is necessary to evaluate integrals of Whittaker functions. This need may 
arise out of the requirement to satisfy normalization or orthogonality conditions. In particular, in the analysis 
of time-dependent Compton scattering, it is necessary to evaluate integrals containing the product of two 
Whittaker W Kt p(ir)-functions, where the variable of integration is the second index fi. This is an unusual 
situation that is not covered by any of the previously known formulas for integrals of products of Whittaker 
functions. The required integrals in the Compton scattering application are members of the general family 

f°° u sinh(27m) P_(l/2 - K - iu) £(1/2 - « + iu) w 

I(S) = / ; 5 W K ,iu{x)W n ,iu[Xo)du, (1) 

Jo s + u z 

where x and xq are real and positive, and s and k are complex. This integral converges for all values of s 
in the complex plane, with the exclusion of the negative real semiaxis, provided that IZen ^ i, |, |,. .., if 
Im k ^ 0. It also converges in the special case s = 0, provided TZen ^ ;>>§>§>••• I n this paper we derive 
several exact formulas for the evaluation of the integral I(s) that fully describe all of the convergent cases. 



II. FUNDAMENTAL EQUATIONS 

We shall begin by briefly reviewing some of the basic properties of the Whittaker functions that will 
be useful in our later work. The Whittaker functions W K ^ L (z) and M KjM (z) are confluent hypergeometric 
functions that are related to the Kummer functions <&(a, b, z) and ^(a, b, z) by 8 ' 9 

M K ^(z) = z» +1 ' 2 e- z ' 2 $(1/2 + M - k, 1 + 2 M ; z) , 

W Kt „(z) = z^ 1 ' 2 e-*' 2 *(l/2 + n - k, 1 + 2/i; z) . ' " ' 
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For small values of \z\, the function M Ktfl (z) is given by the power series 

V(V2 
^ (1 

n=0 

where (a) n denotes the Pochhammer symbol, defined by 9 



M ^ z) - e • 2- (l + 2 M ) n n!' 
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The function W KtfM (z) can be expressed in terms of M Kj/i (z) using 8 
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The integrand in equation Q for I(s) is an even function of u, and therefore we can write 

.sinh(27ra)r(l/2-K-m)r(l/2-K + iu) 
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Next we utilize (JSJ) to express W^^xo) as 



w / \ r(-2i«) , , r(2m) 



which can be rewritten as 

W K> i U (x ) 
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By employing the recurrence formula for the gamma function, z T(z) = T(z + 1), we can obtain the 
form 



W K ,™(x ) = 
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r(l/2 -K-iu) T(l/2 -n + iu) 
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Using this result to substitute for W Kl j u (a;o) m © now yields 
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By utilizing the reflection formula for the gamma function, 

T{l + 2iu)T{-2iu)= , ^ , 
smn(27ru) 

along with the symmetry relation [see Eq. JSJ] 

W K ,i„(a;) = W K .- lu {x) , 

we can rewrite (llUf) as 
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Figure 1. Integration around the closed contour C yields I(s) in the limit r — ► oo [see 
Eq. H16(l ]. In this example, k is a real number in the range 3/2 < k < 5/2, and consequently 
there are two simple poles (P and Q) located inside the contour on the imaginary axis [see 
Eq. I|18p]. The imaginary part of s is less than zero in this instance, and consequently there 
is also a simple pole, R, located at u = i^/s in quadrant I. 

This relation can be split into two identical integrals, and consequently our expression for I(s) can be reduced 
to 



u I\l/2 — hi — iu) 



I(s) = -m — — j — — — — — W K ,- iu (x)M Kt - iu (x )du . (14) 

s + vr 1(1 — 2 iu) 



oo 



III. CONTOUR INTEGRATION 

The fundamental expression for the integral I(s) given by Q is clearly symmetrical with respect to the 
interchange of x and xq. We can use this flexibility to select the arguments of the W and M functions in 
such a way that the integration along the curved portion of the closed contour C in Fig. 1 vanishes in the 
limit r — > oo. By employing asymptotic analysis, we find that this occurs if x max is the argument of the W 
function and a; m i n is the argument of the M function, where 

awn = min(a;, x ) , x max = max(x, x a ) . (15) 

Equation Ijl4(l for I(s) can therefore be recast as the complex contour integral 



where 



I(s) = j L(u)du , (16) 

t , \ ■ u r(l/2 — K — iu) 

L(u) = -ni — ■ — - — — — — — W Kt - iu {x max ) M K iu {x min ) . (17) 

s + u z 1(1 — 2 iu) 

We shall proceed to obtain an exact, closed form expression for J(s) by utilizing the residue theorem to 
evaluate the integral in \i<6\ . 

The integrand L(u) has a simple pole located at u = i\fs, where s/s denotes the principle branch of the 
square root function. This pole is located in quadrant II of the complex u plane if Im s > 0, and otherwise 
it is located in quadrant I. In either case, the pole is contained within the closed integration contour C . 
Additional simple poles are located at the singularities of the function T(l/2 — k — iu), which occur where 
the quantity 1/2 — k — iu is equal to zero or a negative integer. At least one of the poles falls in the upper 
half-plane if TZe k > 1/2. The poles are located at u = u n , where 

K ~\~ r ^j ' n = 0,1,..., [Tie k- i] , (18) 
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and [a] indicates the integer part of a. Note that if IZe n < 1/2, then only the pole at u = iy/s is contained 
within the contour C. 

We can now use the residue theorem to write 

[Were— f] 

I(s) = 2m ^2 Res(w„) + 2tu Rcs(iVs) , (19) 

n=0 

where Res(u*) denotes the residue associated with the simple pole located at it = u*. 

IV. EVALUATION OF THE RESIDUES 

The residue corresponding to the simple pole at u = i^/s is easily computed using the formula 

Res(iVs) = lim (u - iy/s) L(u) , (20) 

which can be immediately evaluated to obtain 

Res(i>/5) = - f r(1 / ( 2 1 ^^ ) W K>V - S (x max ) M K>V - S (x min ) . (21) 

Similarly, the residues associated with the simple poles located at u ~ u n are evaluated using 

Res(u„) = lim (u — u n ) L(u) . (22) 

u — >u n 

Because the poles in this case correspond to the singularities of the function T(l/2 — k — iu), we will require 
evaluation of the quantity 

lim (u — u n ) r [ k — iu] . (23) 

By combining |0J and (|18|) with the recurrence relation zT(z) = T(z + 1), we obtain 
(1 . \ r(l/2 - k - iu + n) iT(l + iu n ~iu) 

J- K — IU \ = — ; 7 = , (24) 

\2 / (1/2 -k — iu) n (iu n - iu - n)„ (u - u n ) 



and therefore 



'1 \ i(-l) n 
lim (u - u n ) T ( - - k - iu j = j — , (25) 



where we have used the fact that (— n) n = (— l)"n! Next we need to evaluate the Whittaker functions 
appearing on the right-hand side of Ijl7|l in the limit u — * u n . Using J5J and (|18|) , we find that 

M K , _ iUn (z) = e- z ' 2 z K - n $(-n, 2k - 2n; z) , 

W K , - iUn (z) = e- z ' 2 z K - n *(-n, 2k - 2n; z) . ' ' 

By employing equations (13.6.9) and (13.6.27) from Abramowitz and Stegun, 9 we can rewrite these 
expressions as 

( a + l)n (27) 
W K ,- iUn (z) = (-l)"n! e~ z ' 2 z^' 2 P^\z) , 

where Pn°\z) denotes the Laguerre polynomial, and 

a = 2k - 2n - 1 = -2iu n . (28) 
Combining lJT7)l. (|2*2*|) . f2*5|) . and l|2*7|). we obtain for the residue 

ResK) = " ! e -(*+*o)/2 {xXo) (^P^ ix) pW(x ) . (29) 

4s — or I (a + n + 1) 
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Utilizing this result along with l|19|l and 1|21[1 . we conclude that 

r(l/2- K-iu)T( 

s + u 2 
2 r(i/2 — K + \/s) 

— 71 r(l + 2^) ^"iV* ( Xmax ) ^K,s/i ( X min) ("3Q^ 

- 4tt 2 e^ 2 ^")/ 2 V — [ 0) P (q) (x) P [a) (xn) 

6 ^ T(a + n + l) As -a 2 " lj " 1 ° j ' 

where a = 2k— 2rt— 1. This previously unknown integral formula is one of the main results of the paper. Note 
that the summation is carried out only if IZe k > 1/2. The integral on the left-hand side of H30|) converges 
for all complex values of s with the exception of the negative real semiaxis, provided that IZe k ^ o j f » §j • • •> 
if Im k / 0. When s = 0, the integral converges provided TZe k is not a positive half-integer. 

A case of special interest can be generated by setting xq — x in (|30[1 . The result obtained is the quadratic 
normalization integral, 

u sinh(27ru) lYl/2 - K - iu) T(l/2 - k + iu) , . s , 
" s + u2 " K, in (*) du 

= tt 2 F(1/2 ~ K ±^f) (a ,) M (x) 

T(l + 2 v / i) K ' VsV ; K ' VsV ; (31) 
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47r e > 

^ T(a + n + 1) As - a 2 

n— 



pM(x) 



which is useful in situations involving the development of a series expansion in terms of a set of normalized 
basis functions. In the following sections, we shall proceed to discuss the limiting behavior of l|30(l observed 
when two of the poles coincide, as well as its relation to formulas appearing in the previous literature. 

V. LIMITING BEHAVIOR 

An interesting situation arises if the quantity 1/2 — k + ^/s is equal to zero or a negative integer, because in 
this case the integral I(s) converges, although the first term on the right-hand side of H30|l formally diverges 
due to the appearance of the factor T(l/2 — k + \fs). This occurs when 

- — ' (32) 



v v 2 

where m is a positive integer or zero. Since y/s denotes the principle branch of the square root function, it 
follows that t/s is located in either quadrants I or IV of the complex s plane, depending on whether Im s 
is positive or negative. Hence Key's > in general, and therefore the function T(l/2 — k + y/s) has no 
singularities unless IZe k > 1/2. The values of m yielding singularities for a given value of k are 

m = 0,1,..., [Ken- -] . (33) 

When s — s m , the divergence of the first term on the right-hand side of (|30|l . containing the factor 
r ( 1 /2 — k + y/s) , is exactly balanced by the divergence of the n = m term in the sum, leaving a finite residual 
quantity. This situation corresponds to a coincidence of the pole located at u — iy/s with the pole located 
at u — u m = i(n — 1/2 — m) [see Eq. (tTH|) ]. In this case the resulting pole has order two. The associated 
residue can be computed by using the standard formula for a second-order pole, but it is more efficient to 
approach the calculation by evaluating equation l|30(l for I(s) in the limit s — > s m . The limiting value of the 
sum of the two divergent terms is given by 

K = lim tt 2 r( ^~ / !^y i) W K ^ (x max ) M KiV5 (x min ) 
1(1 + 2^) 

(* + * )/2 Am! (xxq)(^V 2 



I (A + m + 1) 4s — X z 
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where 

A = 2k — 2m — 1 = 2 y'sm ■ (35) 

Equation l|34|) can be rewritten as 



N 

K = lim - , (36 



where 



TV = 7T 2 (s — S m ) r(1 r / ( 2 1+ 2^ i) (W) M K ,7i (^min) 



(37) 



and 

D = s — s m . (38) 

We can demonstrate that the numerator N vanishes in the limit s — > s m as follows. First we use (0} and 
(I32|l along with the recurrence relation for the gamma function to write 

r p w , A = r(i/2 - k + yi + m) = (yi + ^7) r(i + yi- y^Q 

\2 / {1/2 - k + \fs) m {s - s m ) {\/s - ^/~Sm - rn) m 



1 \ f— l) m 

lim ( a -a m )T[--K + y/e =- — {-2^. (40) 



and therefore [cf. Eq. (J^SJ] 

lim f.s — .s„, 1 F I 

2 

Furthermore, based on l|27|l. 128J1 . and l|35(l . we note that 

(A + lj m (41) 
^ (*) - m! (-l) m e^ 2 z (A+1)/2 P^ A) (z) . 

Taken together, i|37|) . (|40|l . and 141(1 indicate that the numerator N vanishes in the limit s — > s m . The 
denominator Z? also vanishes in this limit, and therefore we can employ L'Hopital's rule to evaluate K by 
writing 

K = lim ^/lim^. (42) 

Since dD/ds = 1 and the second term on the right-hand side of (|37() is independent of s, we obtain 

d 9 , sT(l/2- k + Js) , , 

K = s lim — 7T 2 (S - S m ) r(1 + 2Vi j (w) M KiV ^ (* ml n) • (43) 



Upon differentiation, we obtain after a fairly lengthly calculation 

_ ^ 2 e-(* +3: ">/ 2 (mq)' A+1 )/ 2 mi P£\x) Pr ( n\x ) 

~ r(A + m + l) 

where 7 E « —0.577 is Euler's constant, X — 2k — 2m — 1, 

d 

H = — In [VF K: ^ (x max ) M K: ^ (aw)] 



1 1 m+1 1 

1e + -+H-2tP(\ + 1)-—-+J2 ~ 
A m + 1 ^— ' n 

n—l 



(44) 
(45) 



and 



Mz) = 4- l»r(«) . (46) 



Combining results, we find that in the special case s — s m = (n — m — 1/2) the integral I{s) is given by 
f°° u sinh(27m) T(l/2 - K -iu) T(l/2 - K + iu) . ... 
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where a = 2k — 2n — 1. The allowed range of values for m is given by which indicates that we must 
have IZe k > 1/2 in order for any of these special cases to occur. Note that the singular term with n = m 
is not included in the sum, since that term is contained within K. Equations i|3U|) and 147|) cover all of 
the convergent cases of the fundamental integral I(s). In Sec. VI we present simplified results obtained for 
certain values of the parameters. 

VI. SPECIAL CASES 

The general nature of the expression for J(s) given by Q30|l encompasses many interesting special cases 
involving particular values for the parameters k, s, x, and xq- In this section, we shall briefly discuss a few 
illustrative examples obtained when the first index n is equal to an integer, in which case the general solution 
for J(s) simplifies considerably. For brevity, we shall focus here on situations with s ^ s m . However, we 
emphasize that formulas similar to those discussed below that are applicable to the case s = s m can also be 
obtained in a straightforward manner by starting with (|47|) rather than i|30|) . 

A. k — 

When k = 0, the summation in (|r?C)|l is not performed at all. Making use of the identities 9 

r(--iu\r(- + iu\ = -, — r , (48) 

\2 J \2 J cosh(Tra) ' 1 ; 



and 



we find that (|30[1 reduces to 



sinh(27rw) = 2 sinh(7ru) cosh(7rw) , (49) 



u smliH WQ iu{x)Wo iu{xo)du = | r(V2 + vg Wo v _ {Xmax) M 0tV - s (x min ) . (50) 

This result is convergent for all complex values of s, excluding the negative real semiaxis. Hence the point 
s = is convergent in this case. 

B. k = 1 

When k = 1, there is one simple pole located at sq = 1/4, and we can make use of the identity 

r (-- - iu] r (-- + iu) = — ^ tt- , (51) 

V 2 J \ 2 J cosh(Tru) (1 + 4m 2 ) ' v ' 

along with (|4*9"|) to reduce (|3U1) to the form 

u sinh(7ru) 



2lK ! ^W ltiu {x)W liiu (x )du 

7T r(V5- 1/2) w . , , , . i«„ e -( i+i »)/ 2 



o (1 + 4m 2 )( S + u 2 ) 



- 8 T (l + 2^7s) Wl ^ iXmax) Ml ^ {Xmin) 2 4s -1 

The right-hand side converges for all complex values of s with the exception of the point s = 1/4 [which 
must be treated using l)47p] and the negative real semiaxis. The point s = is convergent. 
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C. k — 2 



In this case there are two simple poles, located at so = 9/4 and si = 1/4. Utilizing the identity 

3 . \ / 3 . \ 16tt 

iu 1 V % 

2 / \ 2 

along with (|49|l , we can simplify i|3U|l to obtain 
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Evaluation of the Laguerre polynomials yields 
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which is convergent for all complex values of s, excluding the negative real semiaxis and the points s = 1/4, 
s = 9/4. These two points must be treated using (|47|l . Note that the point s = is convergent in this case. 
Similar results can be obtained for any positive or negative integer value of k. The integral formula given by 
(|55|l is of particular significance in treating the scattering of radiation in an ionized plasma with a constant 
temperature, as discussed in Sec. VII. 



VII. APPLICATION TO THERMAL COMPTONIZATION 

One of the most important physical applications of the results developed in this paper involves the 
repeated Compton scattering of photons by a hot Maxwellian distribution of electrons with temperature 
T e and number density n e in an ionized plasma. This process, referred to as "thermal Comptonization," 
is the primary mechanism responsible for the production of the radiation spectra observed from celestial 
X-ray sources such as active galaxies, black holes, and neutron stars. 2 When the electron temperature T e is 
constant, the Green's function, / G , describing the temporal evolution of an initially monoenergetic radiation 
distribution satisfies the Kompaneets partial differential equation 1 



dfa 
dy 



_L__d_ 

x 2 dx 



/g+ dx 



where the dimensionless photon energy and the dimensionless time are denoted by 



x(e) = 



y(t) = n e a T c- 



(t - t ) , 



(56) 



(57) 



respectively, and the quantities e, to, t, cr T , m e , c, and k represent the photon energy, the initial time, the 
current time, the Thomson cross section, the electron mass, the speed of light, and Boltzmann's constant, 
respectively. The terms proportional to / G and df a /dx inside the parentheses on the right-hand side of (|56() 
express in turn the effects of electron recoil and stochastic (second-order Fermi) photon energization. At the 
initial time t = t , the radiation distribution is monoenergetic, and the Green's function satisfies the initial 
condition 



f G (x,x ,y) 



y=0 



c 2 5(x 



xo) 



(58) 



where the dimensionless initial energy is given by 

*° = w.- (59) 

By operating on 156(1 with J Q x 2 dx, we can establish that / G has the convenient normalization 

x 2 f G (x, Xq, y) dx = constant = 1 , (60) 



where the final result follows from the initial condition [Eq. I(58|l ] . Note that this normalization is maintained 
for all values of y, which reflects the fact that Compton scattering conserves photons. It can be shown based 
on H56fl that the Laplace transform of the Green's function, 



F(x,x ,s) = I e sy f G {x 1 x ,y)dy , (61) 

is given by 3 



rfc-3/2) 
r(i + 2 M ) 

where the quantity fi is a function of the transform variable s, defined by 



F(x,x ,s) = *o 2 ar 2 e(*°-*)/ 2 ^ , /J M 2 ^(x min ) W 2 ^(x max ) , (62) 



gXl/2 



4) , (63) 
and 

Xmm = mhi(x, Xq) , x max = max(ir, x ) . (64) 
The solution for the Green's function is obtained by performing the inverse Laplace transformation using 
the Mellin integral, 

f G (x,x ,y) = — / e sy F(x,x ,s)ds , (65) 

where the real constant 7 is chosen so that the line IZe s = 7 lies to the right of the singularities in the 
integrand. By transforming the variable of integration from s to 

a' = s + H , (66) 

we can obtain the equivalent expression 

f G (x,x ,y) = —r- / e s yF(x,x ,s')ds' , (67) 

where 

F(x, Xo , s>) = x- 2 x- 2 a(*°-*V 2 3 /_ 2) M % ^ (x min ) W 2>V7 (x max ) . (68) 

1 (1 + lys') 

The exact solution for the Green's function f G (x,Xo,y) can be obtained by taking the inverse Laplace 
transformation of (|55|) . which yields 

1 r +lo ° s „r(Vi-3/2) 
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where we have interchanged the order of integration in the double integral. The inverse Laplace 
transformations on the right-hand side of l|69() are elementary in nature and can be evaluated using the 
formula 

1 p+ico sy 

■ ds = e~ ky . (70) 



27ri ./ 7 _ioo s + fc 



By utilizing this result in l|69[) . we obtain 

1 p+ico r- _ o/oA 

' _ ioo eSV r(i + 2Vi) W ^ s (Xmax) M ^ {Xmin) ds 

32 Z" 00 a u sinh(7ru) ( 71 \ 
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ra e 9y/4 + (2-x)(2-xo)e a / 4 



We can now combine l|67[) . I|68|) . and (|71|l to show that the exact solution for the time-dependent Green's 
function is given by 3 
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Since the fundamental partial differential equation (|56(l is linear, the particular solution for the radiation 
distribution corresponding to an arbitrary initial spectrum can be found via convolution using the Green's 
function. The result given by l|72(l is therefore of central importance in the field of theoretical X-ray 
astronomy. 

VIII. CONCLUSION 

In this paper we have developed several new formulas for the evaluation of a family of integrals containing 
the product of two Whittaker W / Kj ^(x)-functions, when the integration occurs with respect to the second 
index fi, and that index is imaginary. The fundamental integral we have focused on in this paper is 

f°° usmh(2nu)T(l/2- K-iu)T(l/2- K + iu) , . , s , 

I(s)= i '—^ , ? L W K , iu (x)W K , iu {x )du . (73) 



This is related to the Whittaker function index transformation discussed in Refs. 12 and 13. An expression 
of particular interest is the quadratic normalization integral given by l|31|l . The results presented in (|30() 
and H47fl for I(s) allow the exact evaluation of all of the convergent cases of this integral without the need 
to resort to numerical integration. We also point out that by utilizing equations (0, one can easily obtain 
a set of analogous integration formulas applicable to the Kummer functions $(a, b, z) and ^(a, b, z). While 
integrals of this precise type have not been considered before, it is worth noting that I(s) is a member of a 
wider group of integrals containing the product of two Whittaker ^-functions. In general, the other integrals 
in this group involve integration with respect to one of the other parameters, rather than the second index 
as we have considered here. We briefly review a few of these related integrals below. 

Several formulas are available in the previous literature for evaluating the integral of the product of 
two Whittaker T / F K:A ,(a;)-functions with respect to the primary argument x. For example, based upon 
equation (9.12) from Buchholz 10 or equation (20.3.40) from Erdelyi et al. 11 or equation (7.611.3) from 
Gradshteyn and Ryzhik, 8 we have 
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which is valid provided |7£e/x| < 1/2. We note that the formulas in Refs. 8 and 11 are missing a factor 
of tt, and the formula in Ref. 10 contains two incorrect signs. Another closely related example is given by 
equation (7.611.6) from Gradshteyn and Ryzhik 8 or equation (20.3.41) from Ref. 11, 

x^W (x)W (x) dx F(g ± 1} F(g/2 ± 1/2 ± M) r(g/2 ± 1/2 = (75) 
x W R ^{x)W. K ^{x)dx- 2i>/2 + 1 + k)i>/2 + 1-k) ' (?5) 

which is valid provided IZe a > 2 |7£e fj,\ — 1. 
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A few formulas that treat the integration of a product of two Whittaker T4 A K;/J (a;)-functions with respect 
to the first index k have also been known for some time. The most general expression is equation (15.10b) 
from Buchholz, 10 which can be written as 

T(k - iu) T(k + iu) W iUi k _ 1 / 2 (x) W- iu< fc -i/ 2 (^o) du 

/ x ( 76 ) 

= ^ r(2fc) (xx ) k (x + x )- 2fe+1 / 2 K 2k _ 1/2 (~^° * 

where K 2 k~i/2( z ) denotes the modified Bessel function. When k = 1/2, this formula reduces to 
equation (7.691) from Ref. 8, which states that 

[°° S ech(Tru)W iUt0 (x)W- iUt0 (x )du = ^p- e-^)/ 2 . (77) 
Jo x + x 

The new results for 7(s) obtained in this paper [Eqs. (j30|l and l|47(l ] are in some sense complementary to 
these previously known formulas. We emphasize that the expressions developed here are of significance in a 
variety of applications, including the problem of the Comptonization of radiation in an isothermal plasma, 
discussed in Section VII. Our general approach may also allow the determination of the Green's function 
solution for the one-dimensional Schrodinger equation with the Morse potential. 14 We plan to pursue this 
question in future work. 

The author would like to gratefully acknowledge the insightful comments provided by the anonymous 
referee, which led to simplifications in the main derivation and also helped to broaden the applicability of 
the results. 
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